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A component of dark energy has been recently proposed to explain the current acceleration of 
the Universe. Unless some unknown symmetry in Nature prevents or suppresses it, such a field 
may interact with the pressureless component of dark matter, giving rise to the so-called models of 
coupled quintessence. In this paper we propose a new cosmological scenario where radiation and 
baryons are conserved, while the dark energy component is decaying into cold dark matter (CDM). 
The dilution of CDM particles, attenuated with respect to the usual a~ 3 scaling due to the interacting 
process, is characterized by a positive parameter e, whereas the dark energy satisfies the equation 
of state p x = uip x (u> < 0). We carry out a joint statistical analysis involving recent observations 
from type la supernovae, baryon acoustic oscillation peak, and Cosmic Microwave Background shift 
parameter to check the observational viability of the coupled quintessence scenario here proposed. 

PACS numbers: 98.80.Es; 95.35.+d; 98.62.Sb 



I. INTRODUCTION 

According to Einstein's general theory of relativity, the 
dynamic properties of a given spacetime are determined 
by its total energy content. In the cosmological context, 
for instance, this amounts to saying that to understand 
the spacetime structure of the Universe one needs to iden- 
tify the relevant sources of energy and their contributions 
to the total energy momentum tensor. Matter fields (e.g., 
baryonic matter and radiation), are obvious sources of 
energy. Nevertheless, according to current observations, 
two other components, namely, dark matter and dark en- 
ergy, whose origin and nature are completely unknown, 
are governing the late time dynamic properties of the 
Universe. Although fundamental to our understanding 
of the Universe, several important questions involving 
these dark components and their roles in the dynamics 
of the Universe remain unanswered (see, e.g., [l[ for some 
recent reviews). 

Among these questions, the possibility of interaction 
in the dark sector (dark matter-dark energy) , which gave 
origin to the so-called models of coupled quintessence, 
has been largely explored in the literature 0, H| • These 
scenarios are based on the premise that, unless some 
special and unknown symmetry in Nature prevents or 
suppresses a non-minimal coupling between these com- 
ponents (which has not been found - see, e.g., [3] for a 
discussion), such interaction is in principle possible and, 
although no observational piece of evidence has so far 
been unambiguously presented, a weak coupling still be- 
low detection cannot be completely excluded. 
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/From the observational viewpoint, these models are 
capable of explaining the current cosmic acceleration, as 
well as other recent observational results Q. From the 
theoretical point of view, however, critiques to these sce- 
narios do exist and are mainly related to the fact that in 
order to establish a model and study their observational 
and theoretical predictions, one needs first to specify a 
phenomenological coupling between the cosmic compo- 
nents. 

In this concern, an interesting step towards a realistic 
interaction law was given recently by Wang & Meng in 
Ref. [| (see also @) in the context of models with vac- 
uum decay, a class of coupled quintessence in which the 
dark energy equation of state (EoS) is w = —1. Actu- 
ally, in certain sense, one may say that coupled dark en- 
ergy or quintessence models are the natural inheritors of 
the so-called time- varying A(i)-cosmologies @, H, HI J53- 
However, instead of the traditional approach, Refs. [5|, |fj| 
deduced a new interaction law from a simple argument 
about the effect of the dark energy on the cold dark mat- 
ter (CDM) expansion rate. The resulting expression is 
a very general law that has many of the previous phe- 
nomenological approaches as a particular case. 

In this paper, we extend the arguments of Refs. [!, [1] 
to a dark energy/dark matter interaction, where the 
dark energy component is described by an equation of 
state p x = ujp x (w < 0), and explore theoretical and 
observational consequences of a new scenario of cou- 
pled quintessence. Differently from other interacting 
quintessence models, we do not consider interaction be- 
tween the dark sector and the baryonic content of the 
Universe. We also emphasize that this process of inter- 
action is completely different from the physical point of 
view from unification scenarios of the dark sector, an idea 
that has been widely discussed in the recent literature 




We have organized this paper as follows. In Sec. II 
the interaction law and the basic field equations of the 
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FIG. 1: q(z) in the scenario of coupled quintessence, a) Deceleration parameter as a function of redshift for the phantom case 
in which u) = —1.2 and selected values of e. b) The same as in the previous Panel for the quintessence case ui = —0.8. As 
discussed in the text, the effect of a positive e parameter is to decrease the value of q(z), therefore increasing the value of the 
transition redshift zt- q(z) versus redshift for the specific value of e = 0.1 and selected values of EoS parameter. 



model are presented. The influence of the dark energy- 
dark matter coupling on the epoch of cosmic acceleration 
is also discussed. In order to test the observational via- 
bility of the model, Sec. Ill presents a statistical anal- 
ysis involvin g th e most recent type la supernovae (SNe 
la) data fl2l Il3l EH EE EH , observations of the baryon 
acoustic oscillation (BAO) peak (measured from the cor- 
relation function of luminous red galaxies) [t7| and the 
current estimate of the Cosmic Microwave Background 
(CMB) shift parameter from WMAP-5 [H[. In Sec. IV 
we end this paper by summarizing our main results. 



II. THE MODEL 

For a spatially flat, homogeneous and isotropic sce- 
nario driven by matter (baryonic + dark) and radiation 
fields and a negative-pressure dark energy component, 
the Einstein field equations can be written as 



87rG(p 7 + Pb + pdm + Px) = 3-j , 



8ttG(p 7 



Px) 



(1) 



(2) 



where p 7 , pb, pdm and p x , are the energy densities of the 
radiation, baryons, cold dark matter and dark energy, 
respectively, while p 7 and p x are the radiation and dark 
energy pressures. We will always refer to the dark matter 
quantities with the subscript (dm), in order to distinguish 
them to total matter quantities, for which we will use the 
subscript (m). Thus, in our notation, p rn = pb + Pdm- 

By assuming that the radiation and baryonic fluids are 
separately conserved, the energy conservation law for the 



two interacting components (u a T^f 
T& + breads 



= 0, where T a<3 = 



Pdr, 



-J- Pdr, 

a 



~Px ~ %-{Px +Px) ■ (3) 

a 



Now, to complete the description of our interacting 
quintessence scenario we need to specify the interaction 
law. In principle, if the quintessence component is decay- 
ing into CDM particles, the CDM component will dilute 
more slowly compared to its standard (conserved) evo- 
lution, p m oc a~ 3 . Therefore, if the deviation from the 
standard evolution is characterized by a positive constant 
e we may write 1 



Pdm — PdmOd 



-3+e 



(4) 



where e is a constant parameter and we have set the 
present-day value of the cosmological scale factor ao = 1 . 
In what follows we also consider that the dark energy 
component is described by an equation of state p x = ujp x , 
where the constant lu is a negative quantity. 

Now, by integrating Eq. © it is straightforward to 
show that the energy density of the dark energy compo- 
nent is given by 



Px — PxO& 



-3(l+u) 



ePrfmO 3 +e 



(5) 



where the integration constant p x o is the present-day 
fraction of the dark energy density. Clearly, in the ab- 
sence of a coupling with the CDM component, i.e., e = 0, 



1 The positiveness of the interacting parameter e is required from 
thermodynamical arguments. For a discussion, see [6|. 
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FIG. 2: The transition redshift zt in the scenario of coupled quintessence, a) zt as a function of to for some selected values of 
the coupling parameter e. b) The plane (z t — e) for some selected values of uj. In both diagrams, the horizontal lines represent 
the 2a limits on Zt as given in Ref. [lit ]. 



the conventional non-interacting quintessence scenario is 
fully recovered. For oj = — 1 and e =/= 0, we may identify 
p X Q = p V Q (the current value of the vacuum contribu- 
tion), and the above expression reduces to the vacuum 
decaying scenario recently discussed in Refs. [f| [f| . 

Neglecting the radiation contribution, the Friedmann 
equation |T|) for this interacting dark matter-dark energy 
cosmology can be rewritten as 

H = (W 3 + |^« e " 3 + ^a- 3(1+ ^ V2 , (6) 
V 3|w|-e J 

where TL — H(z)/Hq, Qb and fldm are, respectively, 
the normalized Hubble parameter, and the baryons and 
CDM present-day density parameters (for these quanti- 
ties, we have dropped the subscript '0' for convenience). 
The parameter Q x is defined, in terms of the density pa- 
rameter of the dark energy component O x , as 

n^n*--^-, (?) 

6\u}\ — e 

and, therefore 

f2 x = l-0 6 -ff^. (8) 

3\u\ — e 

The above expression clearly shows that the conven- 
tional (non-interacting) quintessence scenario is consid- 
erably modified due to the dark energy decay into CDM 
particles. It is also worth noticing the importance of the 
baryonic contribution to this sort of scenario. Going back 
to high redshifts we see that the presence of an explicit 
baryon term - redshifting as (1 + z) 3 - is well justified 
since the decaying dark enegy component slows down the 



variation of CDM density. Although being subdominant 
at the present stage of cosmic evolution, the baryonic 
content will be dominant (in comparison to CDM) at 
very high redshifts. Actually, it becomes subdominant 
just before nucleosynthesis (z ~ 10 10 for e ~ 0.1), so 
that the CDM component drives the evolution after the 
radiation phase. However, even at late times, baryons 
have several dynamical effects. In particular, they al- 
ter considerably the transition redshift, i.e., the redshift 
where the current accelerating regime begins @ . 

To quantify this latter effect, let us consider the tran- 
sition redshift, Zt, at which the Universe switches from 
deceleration to acceleration or, equivalently, the redshift 
at which the deceleration parameter vanishes. From Eq. 
([T]), it is straightforward to show that the deceleration 
parameter, defined as q(t) = —ad/d 2 , now takes the fol- 
lowing form 

_ 1 1 + 3 -Wf?^ 1 + z r € + (i + -Mfcg + z ^ 
qz ~ 2 i + 3fe%r( 1 + ^ £ + fr( 1 +^ 

(9) 

Some interesting features of the above expression must be 
explored. First, if il x — e — one finds q = 1/2, as ex- 
pected for a flat matter-dominated model (Qb+£ldm = !)• 
Note that the expression of Cl x has been defined by Eq. 
(|8"p . In comparison with the conventional non- interacting 
quintessence scenario, the coupling term (e) modifies 
considerably the transition deceleration/acceleration. In 
principle, since the CDM density scales as a _3+e , while 
the baryon density scales with a -3 , the latter becomes 
dynamically more important in comparison with the non- 
decaying scenario. The overall baryon effect is to delay 
the transition epoch relative to previous cases (including 
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the standard ACDM model), which seems to be in bet- 
ter agreement with some recent SNe results indicating 
z t = 0.46 ± 0.13 at ler Besides this effect, we also 
have in our model the effect of the dark energy EOS oj, 
which introduces a new parameter relative to the decay- 
ing vacuum model @. To better visualize the effect of 
the dark energy EOS, as well as the effect of the dark 
energy-dark matter interaction parameter e, we plot in 
Figure [1] the deceleration parameter as a function of red- 
shift z for some selected values of the EoS (w = —1.1 and 
—0.8) and decaying (e = 0.1) parameters. 

Figures (la) and (lb) show that the effect of a positive 
e parameter, as required by thermodynamical arguments 
p, is to decrease the value of the deceleration parame- 
ter, then increasing the value of the transition redshift. 
The net effect of the dark energy EOS, however, is not 
monotonic as the effect of the interaction. We have in 
general, for values of the dark energy EOS close to the 
vacuum value (u> = —1), that the present value of the de- 
celeration parameter increases for increasing u>, although 
this behavior may change for values of lo much different 
from the standard value, or at high redshifts, as seen in 
Figure (lc). 

Figures (2a) and (2b) show the direct effect of the e and 
u> parameters on the transition redshift (z t ). In agree- 
ment with Figure (1), we see from Panel (2a) that the 
effect of the interacting parameter is in general to in- 
crease Zt- We can also see that the transition redshift 
has a maximum for values of the u parameter a little 
larger than —1. Similar result can also be taken from 
Panel (2b), although z t has a more complicated depen- 
dence on w. We also show in both panels the 2a limits 
given by Riess et al. (2004) [14J . It is clear that values of 
lu larger than —1 as well as smaller than —1 are favoured 
by these limits, whereas larger values of the e parame- 
ter are disfavoured by these limits when u « -1, This 
latter result is in full agreement with the results of our 
statistical analyses discussed in the next section. 



a value in good agreement with the constraints derived 
from primordial nucleosynthesis [l9| . Now, concerning 
the Hubble parameter, it should be recalled that the es- 
timates of Hq through different methods fall on the range 
62-74 km/s/Mpc with an uncertainty of about 10% [2Cl |. 
In what follows, we consider the Hubble Space Telescope 
(HST) Key Project final result, i.e., h = 0.71 ±0.08 
as a Gaussian prior on the Hubble parameter. 



A. SNe la 

The predicted distance modulus for a supernova at red- 
shift z, given a set of parameters s, is 



H P (z\s) = m — M — 51ogG?£ + 25, 



(10) 



where m and M are, respectively, the apparent and ab- 
solute magnitudes, the complete set of parameters is 
s = (H a , tldm, e, w) and stands for the luminosity dis- 
tance (in units of megaparsecs) . 

We estimate the best fit to the set of parameters s by 
using a % 2 statistics, with 



N 

E 

i=l 



[fi P (z\s) -/4(z|s)]' 



(11) 



where ^(zjs) is given by Eq. (fT0|) . /4(z|s) i s the ex- 
tinction corrected distance modulus for a given SNe la 
at Z{, and <Ji is the uncertainty in the individual distance 
moduli. 

In our analysis, we use a combined sample with N = 
192 SNs also used by Davis et al. (2007) [13]. This 
sample consists of the best quality light-curves SNs of 
Wood-Vasey et al. (2007) [l|, which are 60 ESSENCE 
supernovae [13[ |. 57 SNLS supernova e [Tal . and 45 nearby 
supernovae. We also include, as in [12|], 30 new released 
SNe Ia^dassified as "gold" supernovae by Riess et al. 
(2007) [H]. 



III. ANALYSIS AND DISCUSSION 



B. BAO 



The description of the model discussed in the previous 
Section clearly shows that it comprises a multitude of 
cosmological solutions. In a model with such a wealth of 
different possibilities constraints on the parameter space 
arising from current observational data are likely to rule 
out many of the possible scenarios (combinations of e, 
w and £ldm) for the evolution of the Universe. In this 
Section we investigate such observational constraints by 
placing cosmological bounds on the parametric spaces 
e — w and e — Qdm from statistical analyses involving a 
large set of cosmological observations. To this end we 
use the most recent distance measurements to SNe la 
[l2T ] and the current estimates of the baryon acoustic os- 
cillations found in the SDSS data [l7|, as well as, the 
shift parameter from WMAP observations [HI]. In our 
analysis we fix Sit, = 0.0416 also from WMAP results, 



The Baryon Acoustic Oscillations (BAO) given by the 
acoustic oscillations of baryons in the primordial plasma, 
leave a signature on the correlation function of galaxies as 
observed by Eisenstein et al. (2005) [TtJ. This signature 
furnishes a standard rule which can be used to constrain 
the following quantity [l7T ]: 



A 



n 



1/2 



a/3 



-r(z„) 



2/3 



0.469 ±0.017, (12) 



where TL is given by Eq. ((6]), z„ = 0.35 is a typical red- 
shift of the SDSS sample, and T(z„) is the dimensionless 
comoving distance to the redshift z* . As has been shown 
in Ref. [22|, this quantity can be used for models which 
do not have a large contribution of dark energy at early 
times. 
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FIG. 3: The results of our statistical analyses, a) Confidence contours at 68.3%, 95.4% and 99.7% in the plane f2 m - e from a 
joint analysis involving SNe la + BAO + CMB shift parameter + Ho. As discussed in the text, this analysis constrains e to 
values very close to zero (~ 0.09 at 3<r). b) Same as Panel a for the plane ui - e. 



C. CMB shift parameter 



and 



A useful quantity to characterize the position of the 
CMB power spectrum first peak is the shift parameter, 
which is given, for a flat Universe, by (23j : 



ft 



1.71 ±0.03 



(13) 



where z r = 1089 is the recombination redshift and the 
value for ft above is calculated from the MCMC of the 
WMAP 3-yr in the standard flat ACDM model 0]. 

As mentioned above, we also include a Gaussian prior 
on h, as given by the final results of HST Key Project 
[2l| . Thus, in our statistical analysis we minimize the 
following quantity: 



192 

E 



Hth, 



A - 0.469 



0.017 



ft -1.71 
0.03 



ft -0.72 
0.08 



(14) 



0.000 



+0.027+0.057+0.088 
-0.000-0.000-0.000 ' 



D. Results 



with the relative \ 2 l v — 1-03, where v is the number of 
degrees of freedom. These results are much more con- 
straining than those obtained in Ref. [(| for the case 
in = -1, as we have used more recent CMB and SNe la 
data. While in the above reference the bounds on the 
interacting parameter were e = 0.06 ± 0.10 at 95.4% c.L, 
we have found, at the same level, e = 0.000^0 ooo' which 
clearly constrains this parameter to values very close to 
the standard non- interacting case (e = 0). 

In Figure ([3p) we show the plane (u> — e) when the total 
likelihood is marginalized over the density parameter f2 m . 
For this analysis, we have found 

, , _ _i nnK +0. 117+0. 188+0. 258 
UJ — l.UUU_ llg _ .205-0.296 > 

whereas the bounds for e are very similar to those found 
in the previous analysis (Fig. 3a). Clearly, the standard 
ACDM is preferred by this analysis, although much space 
is left for an EoS distinct from — 1. The so-called phan- 
tom models (u < —1) are slightly more favoured by this 
analysis than quintessence (lo > —1) scenarios. 



In Figure (3) we show the main results of our statisti- 
cal analyses. As usual, the total likelihood is written as 
C oc e~ x "/ 2 , where \ 2 is given by Eq. (fl~4"|) . By marginal- 
izing C over the EoS parameter u>, we can quantify how 
much the plane (e-f2 m ) can be constrained by the data. 
The contour levels for this analysis are shown on Figure 
(Hi). At 68.3%, 95.4% and 99.7% c.L, we have found, 
respectively, 

O — n ocq+O. 028+0. 047+0. 066 
"m — U.zoy_ Q26-0. 042-0. 058 



IV. FINAL REMARKS 

The current standard cosmological model, i.e., a flat, 
accelerating Universe composed of ~ 1/3 of matter (bary- 
onic + dark) and ~ 2/3 of a dark energy component in 
the form of the vacuum energy density (A) , is fully consis- 
tent with a variety of observational data. Even so, given 
the complexity of the involved phenomena, it is clear that 
in order to obtain a deeper insight into the nature of the 
dark energy and dark matter, it is worth consider, both 
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from the observational and theoretical viewpoint, more 
complex scenarios as, for instance, models with interac- 
tion between these two components. 

In this paper we have discussed some cosmological con- 
sequences of an alternative mechanism of cosmic accel- 
eration based on a general class of coupled quintessence 
scenarios whose interaction term is deduced from the ef- 
fect of the dark energy on the CDM expansion rate. The 
resulting expressions for the model, parameterized by a 
small positive parameter (e), are very general and have 
many of the previous phenomenological approaches as a 
particular case. In particular, the coupled quintessense 
models proposed here may be thought as a natural ex- 
tension of the decay vacuum scenarios discussed a couple 
of years ago 0, @] . 

By combining the most recent SNe la, BAO, CMB 
shift parameter data and the HST results on Hq we have 
shown that strong constraints can be placed on this kind 
of scenario. We have shown that the free parameters of 
the model are constrained to assume values very close to 
the standard ACDM values, i.e., w ~ — 1 and e ~ 0, al- 



though space is still left for an EoS distinct from — 1 and 
the interacting parameter slightly different from zero. It 
is worth emphasizing that in our analysis the EoS w and 
interacting e parameters have been set as constants. In a 
more realistic case, however, such parameters must vary 
with redshift. The theoretical and observational conse- 
quences of this more realistic interacting w(z) scenario, 
as well as a full comparison with the case discussed in 
the present analysis, will appear in a forthcoming com- 
munication. 
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